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We report elastic and inelastic neutron scattering measurements of the high-Tc ferromagnet 
Mni+^Sb. Measurements were performed on a large, Tq — 434 K, single crystal with interstitial 
Mn content of 5 ~ 0.13. The neutron diffraction results reveal that the interstitial Mn has a mag¬ 
netic moment, and that it is aligned antiparallel to the main Mn moment. We perform density 
functional theory calculations including the interstitial Mn, and find the interstitial to be magnetic 
in agreement with the diffraction data. The inelastic neutron scattering measurements reveal two 
features in the magnetic dynamics: i) a spin-wave-like dispersion emanating from ferromagnetic 
Bragg positions (HK2n), and ii) a broad, non-dispersive signal centered at forbidden Bragg po¬ 
sitions (HK2n-|-l). The inelastic spectrum cannot be modeled by simple linear spin-wave theory 
calculations, and appears to be significantly altered by the presence of the interstitial Mn ions. The 
results show that the influence of the interstitial Mn on the magnetic state in this system is more 
important than previously understood. 

PACS numbers: 75.50.Cc, 75.25.-j, 71.15.Mb, 78.70.Nx 


INTRODUCTION 

Mrii+^^Sb is a high Curie temperature (Tc), highly 
anisotropic, metallic ferromagnet. The observation of a 
ferromagnetic state in Mni+^jSb is unusual in its class 
of materials. Most Mn-alloys are antiferromagnetic, and 
other ‘^d transition metal mono-antimonides are not fer¬ 
romagnetic — CrSb, EeSb and CoSb are antiferromag¬ 
netic, TiSb is paramagnetic and NiSb is diamagnetic [T]. 
The ferromagnetic state in Mni+< 5 Sb is highly sensitive to 
substitutions; for example, Cr-doping quickly tunes the 
system towards an antiferromagnetic state ElE]- In-fact 
substitutions of the cation, the anion, or presence of an 
interstitial can each alter Tc, the magnetic anisotropy, 
or the type of magnetic order BHS]- This potential for 
tuning the properties in Mni+^^Sb, and closely related 
Mni+^^Bi, has attracted considerable attention because 
the materials show promise as alternatives to rare-earth 
containing permanent magnets, and as magneto-optic 
mediums [TOHIH]- 

The spontaneous magnetization in Mni+^^Sb is along 
the crystallographic c-axis at high temperatures. How¬ 
ever, the anisotropy decreases on cooling and passes 
through zero at the spin reorientation temperature, Tsr, 
so that the magnetization is in-plane at low tempera¬ 
tures. In nominally J = 0 poly crystalline MnSb the Tq 
is as high as Tc = 587 K ng. However, single crystal 
studies have been unable to produce the material with¬ 
out the inclusion of interstitial Mn ions [8]. These in- 
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Figure 1. (Color online) (a) Temperature dependence of the 
peak intensity at the (110) (circles) and (10 2) (diamonds) 
Bragg positions across Tc measured on HB-3A. Error bars are 
smaller than the size of the markers, (b) The (10 3) Bragg 
peak measured at 10 K [circles] and 443 K [diamonds] as a 
function of goniometer angle cj, which rotates the crystal in 
the scattering plane, (c) One unit cell showing the low tem¬ 
perature (T < Tsr) magnetic structure of Mni+^Sb. The 
interstitial sites are shown here as fully occupied. 


terstitial Mn ions, Mn2, are present in addition to the 
fully-occupied main site Mn ions, Mnl. The Mn2 ions 
are hosted within the hexagonal NiAs crystal structure 
as shown in Eig. [^c). Compositions in the range of 
S ~ 0.05 — 0.2 have been found to be stable. 

The presence of interstitials significantly alters the 
properties of Mni+< 5 Sb. As S increases, the a lattice 
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parameter and unit cell volume increase, but the c lat¬ 
tice parameter, Tc, Tsr and the total magnetization de¬ 
crease 1310 0. Studies combining magnetization mea¬ 
surements with chemical analysis 0120 ] found that the 
interstitial Mn results in a reduction of the ferromagnetic 
ordering temperature as captured in a simple relationship 
between Tq and 6: 

Tc = (577 - 900 X 5) K. (1) 

The electronic structure and magnetic state of 
Mni+^^Sb have been investigated by a number of authors 
via electronic structure calculations [laillEIHSl]. They 
found that MnSb is described as a metallic system with 
localized magnetic moments, and show that the ferro¬ 
magnetic state is stabilized due to the significant hy¬ 
bridization between the Mn 3d and Sb bp orbitals. All 
these calculations, however, have failed to describe key 
experimental observations such as the behavior of the 
magnetic anisotropy [23] , and the effect of the interstitial 
Mn on the magnetic properties m- Coehoorn et al m 
simply discuss Mn2 as an electron donor, and were un¬ 
able to explain its influence on the magnetic properties. 

Given the importance of the interstitials in tuning the 
properties of Mni+^^Sb, it is surprising that detailed the¬ 
oretical investigations of the effect of Mn2 are lacking. 
Experimental results are also limited, and have focused 
on the empirical determination of the 6 dependencies of 
cell parameters, Tq and Tsr. This is probably due to a 
prevailing view in the literature that there is no magnetic 
moment associated with Mn2. This perspective origi¬ 
nated from the analysis of polarized neutron diffraction 
data [231123], using a model with a complex aspherical 
magnetic form factor for Mnl. A large asymmetry in 
the form factor allows the Mnl site to contribute mag¬ 
netic intensity to Bragg reflections that would otherwise 
be forbidden for Mnl scattering. In these analyses, low- 
Q Bragg peaks had to be excluded to model the data. 
Analysis of powder neutron diffraction data [271 l2H] , x- 
ray Compton profile measurements m, and electronic 
structure calculations m raised questions as to how ro¬ 
bust the aspherical magnetic form factor model is. Ad¬ 
ditionally, a density functional theory (DFT) calculation 
on the influence of Pt doping in MnBi did include Mn at 
the interstitial site, created via a vacancy on a main Mn 
site, and predicted a magnetic moment of order — 3 pe 
associated with the interstitial site m- The absence of 
Mn2 from most calculations may be attributable to the 
inherent difficulty in including a partially occupied, dis¬ 
ordered arrangement of interstitial ions into calculations. 
The role of the interstitial Mn in determining the mag¬ 
netic properties of Mni+^^Sb, therefore, is an open ques¬ 
tion which could provide key insights into the material 
properties and tunability. 

Here we show that the presence of the interstitial 
Mn ion has a strong influence on the magnetic state of 
Mni+^jSb, beyond the role of a simple electron donor. 
We use neutron diffraction to measure a large number 
of Bragg peaks from a ^ = 0.13 single crystal sample at 


temperatures above and below both Tq and Tsr. We find 
that the data is described across the entire Q-range by a 
simple model including a magnetic moment on the Mn2 
site, aligned antiparallel to the moment on the Mnl site. 
We perform DFT calculations including an interstitial 
Mn, and find these to carry a sizable magnetic moment in 
agreement with our diffraction results. The DFT results 
are also consistent with the observed reduction in Tc due 
to the presence of Mn2. We also investigate the magnetic 
dynamics of Mni+^jSb via inelastic neutron scattering 
(INS) measurements performed on a large single crystal. 
We identify a spin-wave-like signal associated with the 
ferromagnetic Mnl ions, and are able to partially repro¬ 
duce its dispersion using a localized-moment Heisenberg 
model. In addition to the spin-wave-like scattering, we 
also identify a broad, intense, magnetic response in the 
inelastic spectrum. This signal is not observed in J = 0 
MnBi, and therefore is likely attributable to the presence 
of the interstitials in Mni.isSb. These results highlight 
the strong influence that the presence of the interstitial 
Mn has on the magnetic properties of Mni+^jSb. 


EXPERIMENTAL TECHNIQUES 

The single crystal of Mni+^jSb used in this investiga¬ 
tion is the same ^ 6cm^ crystal that was studied in 
Ref. [30|. Pieces taken from the large crystal were used 
for diffraction and magnetization measurements. 

To determine the spin reorientation temperature, a 
piece of the crystal was aligned at room temperature to 
within 15° of the easy axis (c-axis) using a permanent 
magnet. The magnetization (M) was then measured in a 
Quantum Design Magnetic Property Measurement Sys¬ 
tem, with applied fields H either parallel or perpendic¬ 
ular to this axis. This measurement clearly reveals a 
spin reorientation temperature of ^ 160 K, as shown in 
Fig-ll a). The same crystal was then used to obtain the 
Curie temperature of 434 K, see Fig. [^b) . This high- 
temperature measurement employed the Sample Space 
Oven from Quantum Design, which utilizes thin quartz 
holders that precluded alignment of the crystal. Here, we 
have defined Tc as the intercept of the greatest tangent 
to the M/H versus T data, with the greatest tangent be¬ 
ing observed at 423 K. According to Eq. Tc = 434 K 
gives an interstitiai content of 16%, consistent with the 
vaiue predicted from the iattice parameters in Ref. [30] . 
With an appiied fieid of 6T, the moment is essentiaiiy 
saturated at 2K, and for H T c we measure a satu¬ 
rated moment of 3.20pB/f-u., where f.u. is a formuia 
unit Mni.isSb. The high-temperature data shows a mi¬ 
nor onset of magnetization magnetization at T ^ 580 K, 
see Fig.j^b). This is iikeiy associated with a smail ciose- 
to-stoichiometric MnSb impurity, or a smaii amount of 
an unknown phase. 

Neutron diffraction measurements were performed on 
the HB-3A four-circle diffractometer at the High Flux 
Isotope Reactor (HFIR) at Oak Ridge National Labo- 
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Figure 2. (Color online) Magnetization data showing the two 
magnetic transitions, with (a) revealing the spin reorientation 
at 160 K by examining data collected on an oriented crys¬ 
tal, and (b) demonstrating the bulk Curie temperature near 
434 K. All of the data were collected while cooling in an ap¬ 
plied field, except for the lOkOe data that were collected on 
the first warming measurement above 300K. 


ratory (ORNL), using a neutron wavelength of 1.003 A. 
Full datasets of more than 190 reflections were col¬ 
lected at temperatures of 10 K (<Tsr), 200 K (Tsr< 
200 K <Tc) and 450 K (>Tc), in a closed-cycle refrig¬ 
erator (CCR). 450 K data were counted for 2 s per point 
and 200 K and 10 K data were counted 5 s per point. In 
addition, a set of 10 peaks, which are predicted by model 
1 (see Results section) to be highly sensitive to the direc¬ 
tion of the interstitial moment with respect to the Mnl 
site moment, were counted for up to 60s per point at all 
temperatures. 


INS measurements were performed on the SEQUOIA 
spectrometer at the Spallation Neutron Source (SNS), 
ORNL with the crystal aligned with (H 0 L) in the scat¬ 
tering plane. Measurements were performed with an in¬ 
cident energy {Ei) of 150 meV with chopper frequency of 
/ = 300 Hz, and with Ei = 60meV and / = 180 Hz. This 
gave energy resolutions at the elastic line of 11 meV and 
4meV, respectively. A CCR was used to reach sample 
temperatures between 10 K and 443 K. Full data sets were 
collected at 10 K, 350 K and 443 K with Ei = 150 meV 
over at least a 130° range of sample rotation angles 
with a step size of 0.5°. Limited angular range data 
sets with 1° step size were collected at 40 K, 100 K, 
125 K, 150 K, 170 K, 200 K, 250 K and 300 K. Measure¬ 
ments were performed with Ei = 60meV at 10 K and 
350 K over an 80° range with a 2° step size. The an¬ 
gular step data were combined, a slight misorientation 
of the crystal was corrected for, and cuts through the 
data were performed using the Horace software pack- 
age [31]. Additional INS measurements were performed 
on the US/Japan Cold Neutron Triple-Axis Spectrome¬ 
ter (CTAX) at HFIR, ORNL, using guide-open-80'-open 
collimation and fixed final neutron energy of 5 meV. 


RESULTS 
Neutron diffraction 

Results from the neutron diffraction measurements are 
shown in Fig.[^ The temperature dependencies of (11 0) 
and (10 2) Bragg reflections across Tq are shown in 
Fig. [^a). These reflections have structure factor con¬ 
tributions from all three atoms in the unit cell. The in¬ 
tensities of these Bragg peaks respond to Tc, indicating 
that they are magnetic Bragg reflections, and confirming 
the identification of Tq from magnetization. Figure [^b) 
shows the (10 3) Bragg reflection at 443 K and 10 K. The 
Mnl site in the NiAs crystal structure has the reflection 
condition I = 2n, where n is an integer. Therefore, the 
(10 3) peak has zero nuclear contribution to its structure 
factor from Mnl, and also zero magnetic contribution if 
the magnetic form factor of Mnl is spherical. Both Mn2 
and Sb ions contribute scattering intensity to the (10 3) 
reflection. 

To determine the origin of the increased intensity at 
low temperature of I = odd, (h — k) ^ 3n Bragg peaks, we 
investigated two possible models for the magnetic state 
in Mni+< 5 Sb for all reflections collected. Both models 
use the NiAs structure, hexagonal space group POs/mmc 
(No. 194), which has two formula units of Mni+^jSb per 
unit cell, with Mnl on the 2a Wyckoff site (0 0 0), Sb on 
the 2c Wyckoff site (| | |), and Mn2 on the 2d Wyckoff 
site (| I |). We use standard Miller indices (h k 1) to 
index the reflections. Note that because of the hexago¬ 
nal symmetry the Miller-Bravais notation can be used to 
illustrate equivalent reflections, given by (h k i 1) where 
i=-(h+k) and permutations of h k i give equivalent re¬ 
flections [32]. Recent x-ray diffraction results from the 
closely related compound MnBi identify a slight distor¬ 
tion from hexagonal symmetry below Tsr, however the 
resolution of our neutron diffraction measurement is not 
sufficient to detect a distortion of this size [33]. The dif¬ 
ference between the two models we investigate concerns 
the magnetic component of scattering below Tq. 

The first model comprises a ferromagnetic arrange¬ 
ment of Mnl moments, with Mn2 moments aligned an- 
tiparallel to Mnl. A depiction of this model is shown 
in Fig. [^c) for T <Tsr, i.e. with spins aligned within 
the a-b plane. For T >Tsr the Mnl and Mn2 moments 
are aligned along the c-axis, but still antiparallel to each 
other. The 3d-bp hybridization is expected to induce a 
small moment on the Sb site antiparallel to Mnl m- 
However, considering the expected moment size and the 
steep magnetic form factor for Sb 5p electrons [2^, the 
overall contribution of Sb to the magnetic diffraction pat¬ 
tern is expected to be small. Therefore we do not include 
magnetic intensity from Sb in the fitting for either model 
1 or model 2, consistent with previous reports [25H?7] . 
As Sb and Mn2 have the same reflection conditions, this 
may have the effect of slightly increasing the moment size 
assigned to the interstitial Mn. 

The model was fit to the data using Rietveld refine- 
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ment in the FullProf software suite m- For a random 
distribution of interstitial Mn ions, the average struc¬ 
ture can be modeled by assuming a uniform distribution 
of Mn2 on every 2d Wyckoff site and scaling the scat¬ 
tering intensity by the occupancy 6. The atomic dis¬ 
placement parameters for Sb and Mn2 were constrained 
to be equal, as they contribute to the same reflections. 
The initial fit was performed against the T = 450 K 
(i.e. T > Tq) data set, allowing us to determine the 
interstitial content of the sample without the influence 
of magnetic Bragg scattering. The standard empirical 
correction for extinction used in FullProf was utilised in 
this refinement m, resulting in an extinction coefficient 
of 580(60) which gives a maximum extinction correction 
of 68% on the strongest Bragg peak. Using this fitting 
procedure the value of the interstitial content was deter¬ 
mined to be 6 = 0.13(1), which is reasonably close to 
the value estimated from comparison of the magnetiza¬ 
tion measurement (see Fig. to Eq. For subsequent 
fits of the magnetic model against 10 K and 200 K data 
sets the interstitial content and the extinction parameter 
were kept fixed. The magnetic form factor of Mn^+ was 
used [26l[35]. For the 10 K data set, the model included 
three domains in equal proportions with spins along a, 6, 
and [-1 -1 0] directions respectively. The results of these 
fits are summarized in Fig. ia), (b) and (c) and Table 
Excellent agreement between data and model 1 is found. 

To highlight the magnetic component of the scatter¬ 
ing, Eig. i e) shows a comparison of the lOK data minus 
the calculated structural component of the scattering, 
with the calculated magnetic contribution to the scatter¬ 
ing. The intensity of the magnetic contribution is seen 
to decrease with Q for both data and calculation, as ex¬ 
pected due to the magnetic form factor. The effects of 
the structure factor and polarization factor, which se¬ 
lects scattering from only moments perpendicular to Q, 
in the neutron scattering cross section are present in this 
data, which is why there is not an exact form factor de¬ 
pendence. In addition, some forbidden Bragg positions 
were measured and appear at zero integrated intensity. 
Attempts to fit an interstitial Mn moment aligned par¬ 
allel to Mnl, instead of antiparallel, were unsuccessful. 
The model with an antiparallel moment on the Mn2 site 
describes the data very well. 

The second model is based on that proposed by Haneda 
et al. [36j for MnAs and discussed in detail for MnSb by 
Yamaguchi et al in which Mn2 does not have a 

magnetic moment associated with it. In this case, they 
explain the observed magnetic scattering at (10 3) and 
other I =odd Bragg positions 12SH2Z1 using a highly as- 
pherical magnetic form factor for the main site Mn. Hav¬ 
ing an aspherical magnetization density breaks the sym¬ 
metry conditions that normally result in systematic ab¬ 
sences for I = odd reflections from the 2a Wyckoff site, 
allowing intensity at I = odd, {h — k) % 3n positions. 
The model includes the combined scattering intensity for 
structural contributions from all ions, plus the magnetic 
intensity from the Mnl ions for all reflections. We have 
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Figure 3. (Color online) Results of model 1 fitted in Full¬ 
Prof compared to (a) 200 K, (b) 443 K and (c) 10 K data, (d) 
Model 2 compared to 10 K data. The data points (red cir¬ 
cles) are the integrated intensities of peaks in scans similar 
to those shown in Fig.[^b) for all measured Bragg positions. 
The results predicted from model 1 (open diamonds), and 
model 2 (crosses) are shown. The difference between the data 
and models (solid lines) have been offset from zero for clar¬ 
ity. (e) and (f) Show the extracted magnetic component of 
the scattering. The data minus the calculated structural part 
of the scattering are plotted (orange circles), along with the 
calculated magnetic component of the scattering for model 1 
(open diamonds) and model 2 (crosses). 


added to model 2 the same extinction correction as used 
in model 1. The overall structure factor resulting from 
the Yamaguchi model, with adjusted interstitial content 
and electron occupancies to match our data, is compared 
to the 10 K data in Fig. [^d). For this model we find 
= 2.49 and = 6.21, compared with = 2.46 
and i?F = 4.55 for model 1 at 10 K (Table [^. The purely 
magnetic component of model 2 is shown in Fig. if), 
along with the difference between the data and calcu¬ 
lated structural component. In this case there is signif¬ 
icant deviation between the model and the data over a 
large range of Q, with the intensities of only some mag¬ 
netic peaks being replicated. Therefore model 2 is not 
as effective as model 1 in describing the results of our 
neutron diffraction experiment. This will be examined 
further in the Discussion Section. 
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Temperature (K) 

a (A) c (A) 

Vol. (A^) 

MmiiI (Fb) 

MMn2 (Fb) 

ext 

scale 


Rf 

450 K 

4.193(3) 5.773(4) 

87.9(1) 

- 

- 

580(60) 346(7) 0.325 3.79 

200 K 

4.195(6) 5.723(8) 

87.2(2) 

2.92(4) [|| c] 

-1.4(3) [|| c] 

580 

368(5) 

2.45 

4.05 

10 K 

4.189(5) 5.690(7) 

86.5(2) 

3.54(2) [L c] 

-2.5(1) [± c] 

580 

372(6) 

2.46 

4.55 


Table L Results from Rietveld refinements fitting model 1 to the data. Vol. is volume of the unit cell, Myini and MMn 2 are 
the magnitudes of the magnetic moments on Mnl and Mn2 sites respectively, ext and scale are the extinction and scale factors 
and and Rf are the statistical agreement factors determined by FullProf. is low for the 450 K data set because of the 
shorter collection time used for this data set. The errors on the moment sizes are the estimated standard deviation calculated 
by FullProf. 


DFT calculations 

To gain a better understanding of the influence of the 
Mn interstitials on Mni+^jSb we performed DFT calcu¬ 
lations including an interstitial in the unit cell (see 133 
for technical details). First we readdress the question of 
whether the interstitial Mn has a magnetic moment, and 
if so how it is directed with respect to the ferromagneti- 
cally aligned moments of the host Mn atoms. To this end 
we considered the 2x2x1 supercell MngSbg depicted in 
Fig. a which contains a single interstitial Mn atom. To 
account for the influence of Coulomb interactions among 
the Mn-Sd electrons we use the PBE-bU approximation, 
and vary the Hubbard U parameter between 0 and 8 eV. 
For all cases, we find that it is not possible to stabi¬ 
lize a configuration in which the interstitial Mn is non¬ 
magnetic. The tables in Fig. show that for U = 0,2 
and 4eV the configuration with Mn2 aligned parallel to 
Mnl is higher in energy than the configuration in which 
it is aligned antiparallel, consistent with the results from 
neutron diffraction. For U = 6 and 8eV, however, the 
configuration with the parallel alignment of the intersti¬ 
tial Mn is energetically favored. It is reasonable to expect 
that the Coulomb d-d interactions are at the lower end 
of the scale, given that the very large wave functions of 
the Sb bp electrons can effectively screen the transition 
metal interactions. Related to this, we find for all cases 
that the system remains metallic, even for a Hubbard U 
of 8 eV. 

Next we examine the effect of the Mn interstitials on 
the ferromagnetic configuration of the host Mn moments. 
To investigate the stability of the ferromagnetic ground 
state (FM) we compare its energy with that of two higher 
energy configurations in which the atomic positions re¬ 
main unchanged, but the host Mn moments assume an 
antiferromagnetic configuration either along the c-axis 
(AFM-c) or the a-axis (AFM-a) as shown in Fig.|^ Com¬ 
paring Fig. I^a) and (c), the presence of the Mn in¬ 
terstitial lowers the energy difference between the FM 
ground state and the high energy AFM configurations. 
The energy difference between the FM and the AFM- 
c configuration reduces quite significantly from 286 to 
212 meV and the energy difference between the FM and 
the AFM-a configuration decreases from 326 to 306 meV. 
The results in Fig. are for a = 2eV, but for a U 
of 8eV we find that this qualitative conclusion remains 
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a(A) 
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Msb(pB) E(meV) 

0 

AP 

4.164 

5.558 

3.208 

-2.570 

-0.105 

0 

P 

4.133 

5.523 

3.104 

0.909 

-0.125 
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2 

AP 

4.280 

5.969 

4.052 

-3.810 

-0.131 

0 


P 

4.311 

5.942 

4.082 

3.915 

-0.123 

281 

4 

AP 

4.338 

6.231 

4.386 

-4.261 

-0.166 

0 


P 

4.352 

6.254 

4.417 

4.284 

-0.155 

18 

6 

AP 

4.388 

6.353 

4.566 

-4.518 

-0.194 

0 


P 

4.393 

6.365 

4.575 

4.483 

-0.199 

-119 

g 

AP 

4.425 

6.443 

4.693 

-4.683 

-0.210 

0 


P 

4.430 

6.441 

4.693 

4.633 

-0.223 

-144 


Figure 4. (Color online) The optimized lattice constants, 
the average host Mnl moment, the interstitial Mn2 moment, 
the average Sb moment, and the total energy per interstitial 
Mn2 atom of the parallel (P) and antiparallel (AP) conhgu- 
ration of Mn2 calculated within the PBE+U approximation 
for U=0,2,4,6 and 8eV. 


unchanged [38]. 

It is important to note that the simulations in Eig.j^a) 
and (c) include the indirect influence of Mn2 via their in¬ 
duced changes in the lattice constants. Erom comparing 
the lattice parameters between the undoped and doped 
system in Eig. ^a) and (c) we see that, just like in the 
experimental observations [9], the Mn interstitials tend 
to increase the a lattice parameter while decreasing the c 
lattice parameter. This leads to the question of whether 
the change in lattice parameters is solely responsible for 
the damaging influence of the Mn interstitials on the EM 
state. To address that question, we present in Eig. |^b) 
the stability of the EM state again, for a system which has 
the strained lattice parameters of the doped system but 
does not include the Mn interstitials. The energy differ¬ 
ences of this undoped strained system shown in Eig.|^b) 
differ significantly from those in the doped system shown 
in Eig. [^c). This illustrates that Mn2 affects the ferro¬ 
magnetic state not just via the change of lattice parame¬ 
ters, suggesting that the magnetic exchange interactions 
between the host and the interstitial Mn moments and/or 
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Without Mn2 

a(A) c(A) MmiiApb) E(meV) per 
Mn2Sb2 

Mn2Sb2 FM 
Mn 2 Sb 2 AFM-c 
Mn4Sb4 AFM-a 

4.189 5.998 4.027 0 

4.189 5.998 3.974 286 

4.189 5.998 3.957 326 


(b) 


Without Mn2 

with strain 

a(A) c(A) MmiiApb) E(meV) per 
Mn2Sb2 

Mn2Sb2 FM 

4.280 5.969 

4.074 

0 

Mn 2 Sb 2 AFM-c 

4.280 5.969 

4.047 

286 

Mn4Sb4 AFM-a 

4.280 5.969 

4.028 

339 



With Mn2 

a(A) c(A) MmiiRm-b) MMn2(M-B) E(meV) per 

Mn2Sb2 

MngSbs FM 

4.280 5.969 

4.052 

-3.810 

0 

MngSbg AFM-c 

4.280 5.969 

4.043 

-3.780 

212 

MngSbs AFM-a 

4.280 5.969 

4.018 

-3.736 

306 


Figure 5. (Color online) The lattice constants, the average of 
the absolute value of the host Mnl moments, the interstitial 
Mn2 moment and the total energy per Mn 2 Sb 2 unit cell of the 
FM, AFM-c and AFM-a configurations for the cases without 
((a) and (b)) and with ((c)) an interstitial Mn2 atom, cal¬ 
culated within the PBE+U approximation with U=2eV. The 
lattice constants and internal parameters for the cases (a) and 
(c) are obtained from optimizing the Mn 2 Sb 2 -FM/Mn 9 Sb 8 - 
FM configuration, respectively. In (b) the lattice constants 
from the configuration with Mn2 are used in the unit cells 
without Mn2, giving effective strained lattice constants. 


doping effects play an important role. Previous attempts 
to explain the magnetic properties of Mni+< 5 Sb just via 
doping effects were not successful m- 


Inelastic neutron scattering 

To investigate the magnetic state in Mni.isSb further, 
we performed INS experiments to probe the magnetic 
dynamics. An overview of the results from SEQUOIA 
is given in Figs, [^a) and which show HL-plane and 
QE'-plane color maps, respectively, as well as constant- 
energy cuts through the data in Fig.j^d). These inelastic 
scattering spectra contain contributions from both mag¬ 
netic and phonon scattering. However, the phonon scat¬ 
tering contribution is effectively negligible in the results 


we present. Mn and Sb have nuclear scattering cross 
sections of opposite sign, therefore their scattering con¬ 
tributions largely cancel-out and as a result the acous¬ 
tic phonon scattering from the sample is weak. Addi¬ 
tionally, the phonon intensity follows a dependence 
and is therefore weak in the low-Q zones we investi¬ 
gate. In the low temperature data the phonon scattering 
is suppressed by the Bose-population factor in the en¬ 
ergy range investigated, whereas at higher temperature 
phonons (attributable to both sample and background 
contributions) simply result in an overall increase in the 
background scattering, particularly at low energies, as 
evident in Fig. [^a). Therefore the results presented are 
attributable to the magnetic dynamics of Mni+^^Sb. 

In Fig.j^a) rings of scattering in the HE plane are seen 
dispersing with energy out of the ferromagnetic positions 
such as (10 0) and (0 0 2); compare to Fig. [^c) which 
depicts the (HOL) reciprocal lattice plane. Only L = 
even positions are magnetic zone centers because the Mnl 
magnetic lattice is half the size of the structural lattice 
along the c-axis, see Fig. [^c). The rings persist above 
both Tsr and Tq. The intensity of the rings is rapidly 
suppressed with increasing Q at all temperatures, and it 
persists up to energies of ^ 70meV (see Fig. Ha)), well 
above the phonon cutoff of the sample ^ 30meV, indi¬ 
cating the magnetic origin of the signal. We observe this 
signal down to low energies in the data from our measure¬ 
ment on CTAX, Fig.[^ This shows that the dispersion of 
the magnetic signal along H can already be observed out 
of (10 0) with an energy transfer of 0.5 meV, see Fig. 
This dispersive signal is reasonably reproduced by a spin- 
wave Heisenberg model for Mnl ions, presented in Fig[^ 
discussed further below, and a qualitatively similar signal 
is observed in MnBi m, which contains no interstitial 
Mn, see Fig. [^b). 

In the inelastic spectra, in addition to the spin-wave 
type signal, a broad signal is evident between the rings 
in the HL plane at all temperatures, see Fig. ga). This 
broad signal was not observed in MnBi, corroborating 
that it is a distinct feature of the excitation spectrum 
of Mni+^jSb, Fig. [^b). The broad signal is centered on 
(0 01) and equivalent positions. (0 01) is not an allowed 
Bragg reflection in this system (see Fig. [^c)), and no 
Bragg peak was observed at this position in any of the 
measurements we performed. Figure [^c) highlights that 
this signal is broad but localized in Q-space, and seem¬ 
ingly non-dispersive. Again, the intensity suppression 
with Q and the energy range of the signal at all temper¬ 
atures both indicate that it is magnetic in origin. We 
observe this signal down to the lowest energies we were 
able to resolve in our experimental set-up on SEQUOIA, 
^ 4meV, but are unable to determine if there is an as¬ 
sociated elastic diffuse signal. This signal was previously 
observed at (101) in the triple-axis-spectroscopy work 
of Radhakrishna and Cable m, but we have now been 
able to map its full Q and E dependence and confirm 
the signal’s magnetic origin, in addition to observing it 
at (0 01). 
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(H 0 0) [r.I.u.] (H 0 0) [r.I.u.] 


(H 0 0) [r.I.u.] 


Figure 6. (Color online) (a) Neutron scattering intensity maps of Mni+^Sb in the HL-plane. The data were acquired on 
SEQUOIA at 10 K, 200 K and 443 K, as indicated, with Ei — 150 meV. Slices were averaged over the energy ranges indicated 
and over ±0.1 reciprocal lattice units [r.I.u.] in the (0.5^ 0) direction (perpendicular to (HOO) and (OOL)). The data on 
interstitial-free MnBi in (b) was taken on the ARCS spectrometer at the SNS, at a temperature of 5K and with Ei — SOmeV, 
courtesy of Ref. [39]. (c) A map of reciprocal space for Mni+^Sb indexed on the structural unit cell. The first Brillouin Zones 
(BZs) of the structural unit cell (dashed line), and the magnetic unit cell of the Mnl ions (solid line) are shown. Structural 
(diamonds and crosses) and Mnl magnetic lattice (crosses) Bragg positions are indicated, with systematically absent reflections 
shown as closed diamonds. 


Further investigation of the scattering centered on 
(0 01) is presented in Fig. 10 Cuts taken along (lOL) 
and (HO 1) directions through datasets collected at 10K, 
350 K and 443 K are shown in Fig. T^a) and (b). As the 
(HOI) direction is the zone boundary for the spin wave 
signal (see Fig. I^c)) the cuts in Fig, [l^b) and the color 
map in Fig. [^cjessentially show the broad signal in iso¬ 
lation. The (1 0 L) direction, however, passes through the 
zone center for both the spin wave and the broad signal, 
therefore both are observed in Figs, p^a) and[^d). We 
show the temperature dependence of the (0 01) signal 
in Fig. [Tq|[c). The integrated intensity was determined 
from cuts made along the (HOI) direction through the 
limited-angular-range data measured at Ei = 150 meV 
for all temperatures. The data were averaged over 10- 
20meV and folded along H and L to improve statistics. 
Two Gaussian functions on a flat background were fit to 
the data, with widths constrained to be equal, but am¬ 
plitudes varying independently, and centers fixed at H=0 
and H=l. The resulting area of the Gaussian centered 
at H=0 was corrected for the Bose population factor, 
[1 — exp(—at each temperature and the re¬ 


sult gives the integrated intensity plotted in Fig. [1 Q][c). 

These results show that the broad, Q = (0 01)- 
centered scattering is influenced by Tsr, Fig. [Tq|[c). 
There is a dramatic increase in intensity of the signal 
on warming past Tsr, which is expected for scattering 
from a transverse magnetic fluctuation as the spins reori¬ 
ent from a5-plane to c-axis alignment, because the mag¬ 
netic cross section depends only on the component of the 
magnetic moment perpendicular to Q. This is further 
supported by the observation that the (0 01) scattering 
is more sensitive to Tsr than the scattering at (±101), 
compare the cuts in Fig. [^b). These observations are 
a further confirmation of the magnetic origin of the sig¬ 
nal, and indicate coupling between this scattering and 
the ferromagnetic state of the system. 


Spin-wave modelling 

We attempted to model the inelastic magnetic response 
of the system using linear spin-wave theory. Isotropic ex¬ 
change interactions including up to sixth nearest neigh- 




















































(H 0 1)[r.l.u.] (1 0 L) [r.I.u.] 

Figure 7. (Color online) (a-c) Neutron scattering intensity 
maps for Mni+^Sb from data taken on SEQUOIA with Ei = 
150 meV at 10 K. The data were averaged over iO.lr.Lu. in 
the two Q-space directions perpendicular to the x-axis in each 
case, (d) Markers show constant-energy cuts through the data 
shown in (a) and (b), averaged over 20-25 meV. Solid lines 
show the result of fitting Gaussians plus a flat background to 
the data. The gray diamonds in the (1 0 L) cut show the region 
of data excluded from the fit, as this region is dominated by 
the broad (0 01) type scattering from the (101) position. 



Figure 8. (Color online) Low energy data measured on CTAX 
along the (H 0 0) direction through ferromagnetic position 
(10 0). Measurements were made with energy transfer of 
0.5 meV at 5, 170 and 300 K, as indicated. The data have 
been normalized to number of counts per minute. The 5K 
data has been multiplied by 5 for clarity. 


bors were the minimum required to model the spiu-wave 
frequeucies. The six uearest ueighbor exchauge parame¬ 
ters, Js, are illustrated iu Fig.|^c). Siugle-iou auisotropy 
was uot iucluded because there was uo evideuce for a 
gap iu the spiu-wave spectrum dowu to at least 0.5 meV 
(<1 % of the baudwidth), see Fig. The iuelastic ueu- 
trou cross sectiou for uudamped spiu waves was calcu¬ 
lated usiug the 1 /S formalism outliued iu Ref. [30] aud 


(a) (b) xio-^ 



(H 0 2) [r.I.u.] (1 0 L) [r.I.u.] 



Figure 9. (Color online) Spin-wave dispersion calculated for 
the H-direction (a) and L-direction (b), with calculated in¬ 
tensities corresponding to (H 0 2) and (1 0 L) directions, 
respectively. The white circles are the points extracted from 
cuts through the data (see Fig.lWd) and text) that were used 
to fit the model parameters. (^ Depiction of exchange con¬ 
stants used in the Heisenberg model, from 1st to 6th nearest 
neighbors for Mnl ions. 


appendix A of Ref. SU For comparison with experimen¬ 
tal intensities, the effects of the magnetic form factor and 
an approximation for the instrumental resolution were 
included in the calculation. We used the magnetic form 
factor for Mn^+ from Ref. [35l The resolution function 
was approximated as a Gaussian in energy with a full 
width at half-maximum of 7.5 meV and the results aver¬ 
aged over the same volume of Q as the experimental data 
(Fig.|7| and averaged over six domains. To determine the 
exchange parameters, the model was fit to the dispersions 
along (10 0), (0 01) and (110) directions, extracted from 
the experimental data via a series of constant-energy cuts 
and fits like those in Fig. [^d). We determined the en¬ 
ergy at the (HO2) zone boundary, 70(4) meV, by fitting 
a Gaussian to a constant-Q cut at (1.5 0 2) from which 
a similar cut at (2 0 2) had been subtracted as a back¬ 
ground. 

The results of this modeling are shown in Fig. [^a) and 
(b), which show reasonable correspondence to the spin- 
wave data presented in Fig. [^a) and (b). However, the 
bandwidth of the dispersion does not appear to be well 
reproduced along the H direction, see Fig.. 0a). Only 
Jq can be uniquely determined from our results. Instead, 
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(b) 


(1 0 L) [r.I.u.] 




200 300 400 

Temperature [K] 


Figure 10. (Color online) (a) and (b) Constant-energy cuts 
averaged over 20-25 meV from Ei = 150 meV SEQUOIA data 
measured at 10K (circles), 200K (triangles) and 443K (di¬ 
amonds). Cuts in (a) were averaged over ±0.1 r.I.u. in per¬ 
pendicular Q-space directions, and successive cuts were offset 
2 X 10“^ for clarity. Cuts in (b) were averaged over ±0.2 r.I.u. 
in perpendicular Q-space directions, (c) shows the tempera¬ 
ture dependence of the signal centered on (0 01), determined 
from (HOI) cuts as described in the text. The plotted inte¬ 
grated intensity has been corrected for the Bose population 
factor. 


the linear combinations 


ex — ± 61/3 

/3 = J4 ± 6 J5 
7 = J2 ± 2 J3 ± 2 J5 


were fit. The values best fitting these measurements 
are aS= 20.0 meV, /3S = - 0.8 meV, jS= 4.8 meV, and 
JqS = 2.9 meV. Fig. shows the calculated frequen¬ 
cies and intensities along two of the measured direc¬ 
tions, for comparison to Fig. [^a) and (b). Previously 
Refs. [32] and [30| reported triple-axis spectroscopy mea¬ 
surements on Mni+< 5 Sb which covered limited ranges of 
Q and not reaching the zone boundary. In Ref [32] 
they fit a Heisenberg model up to fourth nearest neigh¬ 
bors (i.e. J 5 = Je = 0 for comparison to our model) 
to the low energy data measured on a Mni.osSb crys¬ 
tal. They found aS = 8.75 meV, /3S = 2.4 meV and 
76 ' = 5.4 meV. These results extrapolate to a bandwidth 
of 130meV along (HO 2 ), which is much larger than the 
70(4) meV we observe, see Fig. 


DISCUSSION 
Magnetic Structure 

The results of our neutron diffraction experiment on 
Mni.i 3 Sb are well described by the magnetic structure 
shown in Fig.j^c), with ferromagnetic alignment of Mnl 
ions, and Mn2 ions aligned antiparallel to Mnl. There 
appears to be no need to include the model of an aspher- 
ical form factor for Mnl to account for magnetic scatter¬ 
ing at I = odd, (h — k) ^ 3n Bragg peak positions for our 
data set. This is consistent with early neutron powder 
diffraction results on Mui+^^Sb EZl, but in contrast to 
the interpretations of Yamaguchi et al. [23 Reimers 
et al. [26] from their polarized-neutron, single-crystal ex¬ 
periments. 

The single crystals investigated in Refs. [25] and [26] 
were of composition Mni.osSb and Mni.ogSb respectively, 
lower interstitial contents than our crystal, and in both 
cases they excluded some low Q data from their analy¬ 
sis. This may explain the seeming discrepancy between 
their interpretation and our results. In the former case, 
the scattering from Mn2 sites would have been weaker, 
and therefore harder to interpret. The higher Q data 
that they included in their analysis is less sensitive to 
magnetic scattering from the interstitial site, but more 
sensitive to the asymmetric form factor they describe. In 
the case of our data, only model 1 , including an antipar¬ 
allel moment on Mn2, can account for all peaks observed 
below Tc across all Q, see Fig. |^c)-(f). Our data is less 
sensitive to any asymmetry in the Mnl form factor be¬ 
cause it was not a polarized neutron measurement and 
the magnetic contribution at high Q is small. Our results 
do not preclude some asphericity in the magnetic form 
factor of Mnl. It is likely, however, that the aspheric¬ 
ity found in the previous models is too large, as some of 
the scattering intensity should have been accounted for 
by the moment on Mn2. This helps to explain the dis¬ 
crepancies between the Yamaguchi model [23 and results 
from x-ray Compton profile analysis [29] 5 and electronic 
structure calculations m, as well as the failure of the 

model at low Q [2Q1 [23123 • 

The size of the magnetic moments found from our Full- 
Prof refinement of model 1 are consistent with magnetiza¬ 
tion measurements. The saturated moment per formula 
unit was found to be 3.20 //b at 2K from magnetization. 
The results from HB-3A give a total moment per formula 
unit of 3.2(1) /iB at 10 K, taking into account the Mn2 oc¬ 
cupancy of 6 = 0.13(1), see Table [Ij The determination 
of a magnetic moment associated with the Mn2 site is 
further supported by the results from our DFT calcula¬ 
tion. Previous modeling of Mui+^^Sb did not include an 
interstitial Mn, and those that attempted to account for 
it merely adjusted the electron count in stoichiometric 
MnSb models 0 uni EIH211 sni Si- We find a sizable 
moment on Mn2 in the calculation, see Fig. 

Given the presence of a magnetic moment on the in¬ 
terstitial site, in a simple local Heisenberg picture we 
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would expect the antiparallel alignment of the intersti¬ 
tial to enhance the overall ferromagnetic state. However, 
it has clearly been observed that the presence of the in¬ 
terstitial acts to reduce Tq [aiai. Our DFT calculation 
shows that the interstitial Mn degrades the stability of 
the ferromagnetic state beyond the simple influence of 
changing lattice parameters (see Fig. [^, in agreement 
with the experimental observations. This indicates that 
a purely localized moment model may not be sufficient 
to describe the magnetic state in ^ 7 ^ 0 Mni+^jSb. 


Magnetic Dynamics 

Our inelastic neutron scattering measurements of 
Mni+< 5 Sb serve to highlight the influence the intersti¬ 
tial Mn has on the magnetic state. The spin-wave-type 
excitations appear to be less well defined than in the 
non-interstitial containing sister compound MnBi (see 
Fig.§, and are modified such that a simple Heisenberg 
model does not describe the observed dispersion. For 
a local moment system, linear spin-wave theory is ex¬ 
pected to capture the form of the dispersion by inclusion 
of a small number of nearest neighbor exchange param¬ 
eters. Although the observed dispersion (e.g. Fig. [!|(a)) 
appears to be associated with the ferromagnetic wavevec- 
tor, and a similar feature is observed in interstitial-free 
MnBi (Fig. [^b) [39]) , we have shown that the measure¬ 
ment is not fully described by the Heisenberg model for 
localized Mnl spins. Fig. [^ 

For the spin wave calculation along (H 0 2), attempting 
to match both the bandwidth and the gradient of the dis¬ 
persion at lower energies requires a large Jq interaction 
(2.9 meV), which probably reflects the inadequacy of the 
model. The shape at the zone boundary remains quite 
different, resulting in the 70(4) meV bandwidth not be¬ 
ing reproduced, see Fig. [^a). The calculation does not 
produce a signal emanating from (0 0 1 ) as is observed in 
the data, see Fig. [^c). As (HOI) is a zone boundary 
for the Mnl magnetic lattice, as illustrated in Fig. [^c), 
the spin-wave intensity at (0 01 ) is concentrated around 
the zone boundary energy E = S5meV (see Fig. |^b)). 
Our crude attempts to extend the model to include an 
interstitial Mn2 moment failed to better reproduce the 
spin wave data, nor did they produce a signal like the 
observed scattering along (0 01). These attempts are 
limited at the present time by the absence of general the¬ 
oretical and computational tools to include the effects of 
a disordered component in a system on the inelastic mag¬ 
netic response. Including a small moment at one inter¬ 
stitial site in our model calculation artificially introduces 
a symmetry element which is not present in the real sys¬ 
tem. Extending the calculation to include thousands of 
unit cells, which can then be randomly assigned inter¬ 
stitials, is currently difficult to implement and extremely 
computationally intensive. 

In addition to the spin wave, we identified a second 
magnetic signal in the inelastic spectrum which is broad 


in Q, and centered on an unexpected position in recip¬ 
rocal space Q = (0 0 1) (see Figs. ilzl and [1Q| ) . Like the 
spin-wave scattering, this signal responds to Tsr, but un¬ 
like the spin-wave scattering, no similar signal is observed 
in interstitial-free MnBi. (0 01) is not an allowed nu¬ 
clear or magnetic reflection (for either magnetic structure 
model discussed above). The presence of a (0 01) Bragg 
peak would require the correspondence between the up¬ 
per and lower halves of the unit cell shown in Fig. [^c) 
to be broken, for example, if Mnl was antiferromagnet- 
ically aligned along c, or if Mn2 was only present in the 
lower half of the cell. We do not observe a Bragg peak at 
(0 01 ) in any of the neutron data sets collected, implying 
that no such long-range order is present. However, the 
periodicity of the signal in reciprocal space implies that 
underlying long-range correlations are present. 


Our results show that the dynamic correlations ob¬ 
served at (0 01 ) are not associated with a long-range 
order in the system. They could, however, be associ¬ 
ated with a short-range structural or magnetic order in 
the system. This proposition may be supported by the 
broad appearance of the signal in reciprocal space (see 
Fig.§, which suggests a short correlation length in real 
space. We made a rough estimate for the correlation 
length by fitting Gaussians to cuts similar to those in 
Fig, [l^a) and (b), integrated over 10-15 meV, from the 
Ei = 60meV, T = lOK data. We find correlation lengths 
of —16 A in the a-b plane, and in the range 10-18 A along 
c (where it is harder to define, due to the overlap of the 
signal with the spin-wave-type scattering), i.e. -4 unit 
cells in plane, and -3 units cells along c. We were unable 
to determine whether there is an elastic, diffuse scatter¬ 
ing signal at the (0 01 ) position due to the neutron in¬ 
strumentation used in this investigation. We do observe 
the signal down to the lowest energies that were probed 
on SEQUOIA, ^ 4meV. A short-range order could be 
present in the system if the interstitial Mn ions are not 
randomly distributed, but instead are structurally or¬ 
dered over a few unit cells. An observation of elastic, 
diffuse scattering at (0 0 1 ) would establish that short- 
range order is present. Alternatively, the presence of in¬ 
terstitials could modify the interactions of neighboring 
main-site Mnl ions, with the associated magnetic excita¬ 
tions unable to propagate over large length scales due to 
disruption from a random distribution of magnetic Mn2 
ions. Eurther studies, including diffuse neutron scatter¬ 
ing experiments, are highly desirable to investigate the 
potential short-range order in this system. 

The spectral weight of the (0 01) scattering, i.e. the to¬ 
tal intensity of the signal integrated over QE space, can 
in principal be used to indicate the strength of the mag¬ 
netic excitations and therefore the size of the fluctuating 
magnetic moment, which could give an indication of the 
origin of the scattering. Unfortunately, due to the overlap 
of the two magnetic excitations in QE space, we cannot 
quantitatively compare the spectral weights of the two 
signals. Qualitatively, however, by inspecting a series of 
slices and cuts through the data, such as those shown in 
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Figs. HI [3 and we see that the strength of the signal 
emanating from (0 01) is comparable to the spin-wave 
signal. This implies that the (0 01) signal is associated 
with a relatively large magnetic moment per formula unit 
in the sample, of the order of the Mnl moment size. This 
makes it unlikely to be purely the result of the Mn2 mo¬ 
ments behaving independently of Mnl. Identification of 
the origin of this signal would be extremely significant 
for understanding the magnetic state in Mni+^^Sb and 
determining the possibilities of tuning the state with in¬ 
terstitial ions of Mn or other 3d transition metals. 


CONCLUSIONS 

We have shown that a magnetic moment on the inter¬ 
stitial Mn in Mni+^jSb is required to describe the neu¬ 
tron diffraction data from a single crystal of Mni.isSb. 
This magnetic moment is aligned antiparallel to the main 
site Mn magnetic moment. We performed DFT calcu¬ 
lations which find that the interstitial Mn is magnetic, 
consistent with the diffraction data. The DFT results 
are also consistent with previous experimental evidence 
that the interstitial reduces Tq and c lattice parameter, 
and increases a in Mni+^jSb. We find that the presence 
of the magnetic interstitial Mn also has a substantial ef¬ 
fect on the magnetic excitation spectrum of Mni+< 5 Sb, 
measured by inelastic neutron scattering. It results in 
the appearance of an intense, broad signal in Q-space, 
that is likely associated with short-range order. While a 
Heisenberg Hamiltonian model calculation captures the 
essential Q-space dependence of the spin-wave type signal 
observed, the additional broad signal cannot be explained 
by straightforward attempts to include an interstitial ion 


in this model. 
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